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Collisionless damping of a condensate of cold Fermi atoms, whose scattering is controlled by a 
Feshbach resonance, is explored throughout the BCS and BEC regimes when small perturbations 
on its phase and amplitude modes are turned on to drive the system slightly out of equilibrium. 
Using a one-loop effective action, we first recreate the known result that for a broad resonance the 
amplitude of the condensate decays as t -1 / 2 at late times in the BCS regime whereas it decays as 
t _3//2 in the BEC regime. We then examine the case of an idealized narrow resonance, and find that 
this collective mode decays as i~ 3//2 throughout both the BCS and BEC regimes. Although this 
seems to contradict earlier results that damping is identical for both broad and narrow resonances, 
the breakdown of the narrow resonance limit restores this universal behaviour. More measureably, 
the phase perturbation may give a shift on the saturated value to which the collective amplitude 
mode decays, which vanishes only in the deep BCS regime when the phase and amplitude modes 
are decoupled. 



I. INTRODUCTION 

We now have remarkable experimental control over 
cold alkali atoms interacting through a Feshbach reso- 
nance. Manipulation of the binding energy through ex- 
ternal magnetic fields enables us to evolve them con- 
tinuously from the weakly coupled BCS-like behavior of 
Cooper pairs to the strongly coupled Bose-Einstein Con- 
densation (BEC) of molecules [I]. The transition is char- 
acterised by a crossover in which, most simply, the s-wave 
scattering length as diverges as it changes sign [5] . 

Recently, a considerable theoretical effort has been ex- 
pended on understanding how such a tunable supercon- 
ductor/condensate responds to an initial non-equilibrium 
change in its order parameter. The behavior of tun- 
able gases can differ according as the Feshbach reso- 
nance is broad or narrow, the former essentially describ- 
ing a one-channel system, the latter a two-channel sys- 
tem. Nonetheless, in [4] it is argued that, for linear per- 
turbations, the relaxation through damping of this col- 
lective motion is the same for both broad and narrow 
resonances. What determines the nature of the damp- 
ing is whether the chemical potential is positive (BCS 
regime) or negative (BEC regime). In the former case it 
is proposed that the perturbed amplitude of the conden- 
sate decays with elapsed time t as t -1 / 2 , as happens for 
superconductors [5], and in the latter as t~ 3 ^ 2 . Advanced 
detection techniques certainly permit experimental time- 
resolved studies of the collective pairing mode [3] with 
regard to such perturbations |4j. [There has also been 
work on non-linear perturbations [6-8], but this goes 
beyond the analysis that we shall present here.] 

Here we provide an extensive study to reexamine the 
relaxational dynamics of the amplitude mode of the 
fermionic condensate as a result of this (Landau) damp- 



ing throughout both BCS and BEC regimes for both 
broad and narrow resonances, building on our earlier 
work [51 [TU] on the BEC/BCS crossover in the pres- 
ence of a Feshbach resonance. In [9] [10] we have shown 
that many condensate properties at T = can be de- 
rived easily from the derivative expansion of the renor- 
malized one-loop effective action. In particular, under 
long wavelength and small energy approximations this 
action can be used to obtain a hydrodynamic description 
of the BEC-BCS crossover from which the equation of 
state, intimately related to the speed of sound, can be 
derived. In this paper we particularly build on the work 
of [5] , creating an extended full one- loop effective action 
to explore the real-time evolution of the collective mode. 



We find that the situation is more subtle than as pre- 
sented in [4 , with some differences in the relaxational 
dynamics of the amplitude mode for broad and narrow 
resonances in the BCS regime with positive chemical po- 
tential. For broad resonances we recreate the results of [4] 
directly. However, for idealised narrow resonances we 
find a f -3 / 2 decay across both BCS and BEC regimes. 
Nonetheless, we anticipate that for more realistic reso- 
nances the long-time behaviour of [4] will be recovered 
in the deep BCS regime, albeit with i -3 / 2 transients. A 
new result is that, even under linear perturbations, the 
amplitude mode is found to decay to a shifted value from 
its initial condition, due to the phase-amplitude coupling, 
irrespective of whether the resonance is broad or narrow. 
Such a shift is expected for non-linear perturbations but 
not for linear perturbations, when typically the system 
relaxes to its initial equilibrium state. These findings can 
be tested experimentally. 
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II. EFFECTIVE ACTIONS 

We consider a condensate comprising a mixture of 
fcrmionic atoms and molecular bosons, in which the 
fermions ip a (x), with spin a = (t,|), undergo self- 
interaction through an s-wave BCS-type term, and two 
fermions can be bound into a molecular boson <fi(x) 
through a Feshbach resonance. The general 'two-channel' 
microscopic action is given by (U > 0, g fixed) 



5' 



J L T.4- 



i d t 



2m 



/' 



+ U i%(x) i/>l(x) i/>±(x) ip t (x) 
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Thus, the bound bosons of Feshbach resonance (field <p) 
have twice the mass of the fermions, M = 2m, and a 
tunable binding energy, v. 

On introducing the auxiliary field A{x) = 
Uil)^(x)'ip^{x), a Hubbard-Stratonovich transforma- 
tion leads to an effective Lagrangian density S quadratic 
in the Fermi fields. We then integrate them out [S] to 
write S in the non-local form 



Snl = —iTrhxQ' 



,11 ,/x<! - ^|A| 2 



V 2 

ldt+ 2M +2 »- V 



*)>, (2) 



in which Q 1 is the inverse Nambu Green function, 

id t — e A(x) 
A* (a;) id t +e 



where 



A(x) = A(x) - gcj)(x) 



(3) 



(4) 



that represents the two-component combined condensate 
(and £ = — V 2 /2m — fj,). 

The combined condensate amplitude and phase of 
A(.t) = |A(x)| e l0 ^ x ^ can be determined from those 
of A(ar) and c/>(x), defined respectively by A(x) — 
\A(x)\ e i6 ^ and <j){x) = -\<j>{x)\ e t6 ^ x \ The U(l) 
invariance of the action under 6* a — > 6 A + const., 6$ — > 
00+const. is spontaneously broken. Thus, 5 Sjvl = per- 
mits spacetime constant gap solutions |A(x)| = |Ao| 
and \<b(x)\ = \4> \ £ (whereby (|A(x)| = |A | + 0). 

We now consider the fluctuations around the gap con- 
figurations and simultaneously perturb in the derivatives 
of 6a and 9$, small perturbations in the scalar conden- 
sate densities QI] 5|A| = |A| - |A | and 6\<f>\ = \<p\ - \<f> \ 
and their derivatives. To guarantee Galilean invariance 
we expand [SI [H]in £ = Q^ 1 — G^ 1 , where Q^ 1 is ob- 
tained from Q~ x by substituting A(x) with the Galilean 



scalar |A(x)|. We construct the condensate effective ac- 
tion S e g at second order in S. Diagrammatically, this 
amounts to taking account of fermionic one-loop effects 
in the effective action. 

S e s comprises the integral of a local density together 
with non-local fermionic cut contributions. Later we 
shall restrict ourselves to homogeneous collective modes. 
For their dynamics we only need the quadratic part of the 
local density, while retaining the full cut contributions. 

(2) 

The relevant S^ s then takes the form 



S^=J^-^\Y6^- ~-^(Y6.,) 

1 



m 



^{6^-6^-2^8^ 64, 
+(2 M - v^m? + jfS\A\6\<f>\ - ^(5\A\) 2 



I d 4 q I d 4 P |Xe AeA (g,p)e A (q)e A (- (? ) 

+M AA (p,q)~6\A\(q)~6\A\(-q) 

+M 0xA (p,q)e A (q)5\A\(-q)\, (5) 



for VW-functions to be given later. In (|5|) <5|A| = 
|A| — |Ao|, with Fourier transform <5|A| and A is 
the Fourier transform of the Galilean invariant Q A = 
6 A + (\76 A ) 2 /4m. The varying coupling strength (as 
v varies with external magnetic field) is U e f t = U + 
g 2 /(2/i — v). The phases of A and <f> axe coupled with 
strength O 2 = (2g/U)\4> Q \\A Q \. Defining E 2 = e 2 + |A | 2 
and e p — p 2 /2m — /1, the fermion number density is 
Po = Po + Po, where p£ = J d 3 p/{2n) 3 [1 - e p /E p ] is 
the explicit fermion density, and Pq = 2\tfio\ 2 is due to 
molecules (two fermions per molecule). For more details 
see 0. 

The separation of the action ^ into local and seem- 
ingly non-local parts is somewhat misleading. The Ms 
have g-independent terms that lead to further contribu- 
tions to the local part of the action (the effective local 
Lagrangian). These are proportional to 9. and (<5|A|) 2 
and have been displayed elsewhere [S] . The remainders of 
the double integrals comprise different non-local combi- 
nations of two-Fermion cuts. They are the source of the 
Landau damping, which occurs in a collisionless regime 
via direct dissipationless energy transfer from the collec- 
tive mode to single particles. For the relatively simple 
.Ms for the problem in hand of a homogeneous system 
it is easier not to make this separation into local and 
non-local parts. 

The dynamical equations of the collective modes can 
be obtained by taking the variation of the above effec- 
tive action with respect to the associated field variables. 
In what follows, we will consider broad and narrow reso- 
nances in turn. 
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III. BROAD RESONANCE 

In the case of a broad resonance the one-channel model, 
which is obtained by eliminating the molecular bosons 
from the above two-channel effective action, suffices to 
describe its dynamics |I2[ I13| . Thus, the corresponding 
action is given by setting <fi = in the action which, 
in turn, gives A — A. 

Consider a spatially homogeneous condensate. If, at 
time t = we turn on a small instantaneous homogeneous 
change in the external field this will induce a homoge- 
neous perturbation <5a in the gap parameter (condensate 
amplitude). On the other hand, since the gradient of 
the phase is the fluid velocity, a spontaneous translation 
(kick) of the system will induce a perturbation 8e A in 
the phase, which we need to include for completeness. 
Although we have not included a trap to constrain the 
condensate in the analysis, to keep the results as sim- 
ple as possible, this translation of the system could be 
implemented through a kick on the trap. 

The linearized equations of motion are Fourier ana- 
lyzed in terms of frequency u) or, more conveniently, ex- 
pressed in terms of the Laplace variable s = — iu> as 



s 2 Me A e A (s 2 ) 6 A + S M 9a a(s 2 ) S\A\ 
-sM e&A {s 2 ) A + Maa(s 2 ) 5\A\ 



$e A 
5a 



(6) 



(where we have dropped the tildes that normally repre- 
sent transforms) . After some straightforward but tedious 
calculation the Ms are found to be 
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The solutions of the Laplace transform of the equations 
are 



Oa(s) 



-Maa(s 2 ) S 9a + sM 8& a(s 2 ) S a 
s 2 V{s 2 ) 



$|A|(a) = s2m ^Sa(s 2 )Sa - sM e& A(s 2 )Se A ^ 



s 2 V(s 2 ) 



where 



V(s 2 ) = M b ^ & (s 2 )Maa(s 2 ) - M 2 e&A (s 2 



(9) 



The real-time evolution of 9(t) and <5|A|(i) can be ob- 
tained from carrying out the inverse Laplace transform 
along the contour in the s-plane shown in Fig.l [14] . 

The singularities include the pole of the phonon mode 
at s — iO as well as the branch cuts, extending from 
s = ±iE t h to s = ±ioo due to the possibility of breaking 
a condensate into fermionic excitations. The solutions to 
equations pi are obtained as 
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V(0) dA 7T 
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A0 A R 



e A e A r 



2?>( W )+2??H 
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dwV R (u)Me A Ai{u) -V i (w)M 6a ar{u) . , . 1JC 

r>2 i \ t,2( \ sm[ut\de A 

B u V 2 R {uj)+Vj{uj) 



(10) 



(11) 



The first terms of the above results come from the contri- 
bution of the pole and the remaining terms are due to the 
contributions of the cuts. All of the Ms and V acquire 
an imaginary part as s = iu crosses the cuts, where their 
real and imaginary parts are denoted by 



M{s = iuj^e) = Mr{uj) ± iMi(ui) , 
V{s = iLuTt) = V R (u)±iD I (u>) . (12) 



From Eq. (JoJ) Vrj are found to be 

V R {u) = [M 6a 6 a rMaar - Me A e A iM A Ai 
-M 2 9aAR +M 2 0&ai ] (w) 

X>j(w) = [Me A e A rM aa i + Me A e A iMAAR 

-2Me A ARM9 A Ai]H- (13) 

Upon turning on the perturbation Sg (Sa), since the 
phase mode is generally coupled to the amplitude mode, 
it can drive not only the phase (amplitude) modes but 
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FIG. 1: This Fig. shows the analytic structure of the two- 
Fermion amplitudes discussed in the text in the complex-s 
plane. The contour is deformed from the Bromwich contour 
that encompasses all singularities in a counterclockwise way. 
In addition to the pole at s = iO due to the phonon, the cuts 
are displayed by a wiggly line. 



also the amplitude (phase) modes away from their ini- 
tial equilibrium values. Because of this coupling the per- 
turbed modes relax to their respective saturated values 
given above. However, since the phase fluctuations can 
couple to the density fluctuations, we can observe this 



collective phase mode from the spectrum of the density- 
density correlation function (see Ref. [15] for details). 
Here we merely focus on the real-time behavior of the 
collective amplitude mode that can be observed by exper- 
imentally time- resolved techniques [3J. The main results 
below follow directly from Eqs.(ll I. 



A. The BCS regime 

In the BCS regime for ^ > 0, the threshold energy is 
E t h = 2|A |. The real and imaginary parts of the Ms 
are summarized as follows: 
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We define w = u)/ep,fj, = fi/ep, |Ao| = |Ao|/cf, and the 
symbol "P means that the principal value of the integral 
is taken. The damping dynamics of the collective mo- 
tion, particularly at late times, largely depends on the 
behavior of the M(uj)s for oj near the threshold energy, 
i.e. uj — E t h + + . These are found to be 
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with 



Vi = 16Me A e A RMg A g A i - 2Me A ARMg A Ai (20) 
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As expected, A^e A Afl and M.e A AR. vanish in the deep 
BCS regime in which fi 3> |Ao|, resulting from particle- 
hole symmetry. The well-known divergence on the BCS 
density of the state at threshold energy E t h = 2|Aq| [12] 
renders Mg A g A i{uj = E t h + + ) singular. Together 
with the behavior of A^aa _/?(/) ( w = E t h + + ) given by 
Eqs.(17), which shows their vanishing at threshold en- 



ergy, it gives 2?r,j, defined in the expressions (13), as 
follows: 



V R (u = Et h + + ) ^V R 



V I (u J = E th + 0+)~V I [ — -1 

^th 
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where 



V R = -8M 2 eA g Al ~M 2 g AAR ; 



As seen from Eq.Q, the pole at s = iO seems fictitious 
with respect to the amplitude perturbation [4], only be- 
coming a true pole when an additional phase perturba- 
tion is turned on. It is the terms involving the singu- 
lar behaviour of Aig A g A j at cj near the threshold energy 
that give the dominant contributions to the late-time dy- 
namics of the amplitude mode. Thus, one finds that the 



integrands in Eq.(ll) behave as 



'E , RMg A 8 A i — T>jMg A e A r 
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f> R Mg A AI-VlMg A AR ijfa) 



These two terms in turn will lead to the different damping 
behaviour after carrying out the integration over ui. 

Putting all this together, we find the late-time solu- 
tions as 
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In the BCS regime, <5|A|(t) decays at late times domi- 
nantly as i -1 / 2 [H[S] due to the amplitude perturbation 
while it has a subdominant decay as i -3 / 2 given by the 
perturbation 6g. 

Additionally, if it were possible to implement small 
phase perturbations, <5|A|(t) decays to a non-zero value 
proportional to the perturbation 5g, resulting from the 
amplitude-phase coupling. This looks to give a shift on 
the saturated value that the amplitude mode decays to 
from its initial value, but this becomes vanishingly small 
in the deep BCS regime where Mg A A — 0. We might 
hope that this subdominant behaviour could be made 
visible by judicious choice of initial perturbations, but 
we do not know. 



B. The BEC regime 

In the BEC regime the chemical potential changes its 
sign, fi < 0. The threshold energy now becomes E t h = 
y // 2(\/i\ 2 + |A | 2 ). The real and imaginary part of the 
Ais obtained from Eqs.(12| are given by 



Me AdAR {uj) = ^(2™e F )5|Ao| 2 / 1 (w,-|/j|,|A |); 

|A | 
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4|o,|yW2) 2 -|A | 2 
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together with Eqs.( 14 ),( 16 ). For lj near the threshold 



energy, the .Ms can be simplified as 
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which, in the deep BEC regime, becomes 

,(2-,/2>|Ao| 2 



M 6a ai - 



m .„ , 1 7T |An| 2 

'2^ (2 ^ )3 8 Vf ; 

7T |Ao| 



Z7T Z 



2V2 H 
ttIAqI 



-2^ (2m ^ 4 |^ 



(27) 



Although the behaviour of the above .Ms seems different 



from that in the BCS regime, "Dfirn follows Eq.( 19 1 with 
^Dr(i) given respectively by 



V ' = 8 r^ M e A e A RMe A e A i - 2M 0& ARMe A A[^ 



So, as lj is near the threshold energy, the integrands in 
Eq.(ll) can be approximated by 
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77Z 



The same damping behaviour will be seen from the above 
two respective contributions. 

Putting all these results together now gives the late- 
time result: 
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Note that E th = \/2(|/i| 2 + |A | 2 ) in the BEC regime. 
The amplitude mode decays as t~ 3 / 2 in both oscillating 
terms, but we note that the coefficient of Sg is the larger 



in the deep BEC regime, as seen from Eqs.(27) in which 
Me A A >• M.e A e A for \fx\ ^> |A |. Insofar as it is possible 
to implement an initial phase shift <5|A|(i) then decays 
to a saturated value with a nonzero shift from the ini- 
tial condition. The significance of introducing the phase 



perturbation is seen from the two results above. 

In summary, our calculation of the generic i -1 / 2 damp- 
ing in the BCS superconductor regime and the t -3 / 2 
damping in the BEC molecular regime reproduces the 
results of [4 obtained by independent methods. As in 
j3], we have seen that the transition between the two 
regimes of the different damping behaviours occurs at 
fx = 0. However, for 5g A ^ 0, there is a residual displace- 
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ment of <5|A|(f) in the BEC regime. 



IV. NARROW RESONANCE 

The situation is rather different when the Feshbach 
resonance of bound bosons is very narrow. We now have 
a two-channel model in which the resonance has to be 
taken into account explicitly. To illustrate this, for sim- 
plicity we consider the idealised case [16] of U = in 
the above Lagrangian where the effects from the back- 
ground scattering length a^ g are ignored as long as the 
system is away from the deep BCS regime [T7]. It leads 
to A = — gcf> from Eq.Q, so the field variables of the col- 
lective modes then are the phase 9<p and the amplitude 

On perturbing the system the collective phenomena 
(in the same notation as before) are described by the 
following equations: 



M e<t ,4,R(uj = E th + 0+) ~ Mg^R 
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whereas the behaviour of the other .Ms remains the same 



as in Eqs.(17) and (18) in which the labels 9 a and A 
are replaced respectively by the corresponding 9g and 
4>, and |A | is changed to g|0ol- F° r 2/i ^ v we have 
■M</></> r(u = E t h + + ) 7^ 0. This is crucial and leads to 
the rather different results for T>rj, which are obtained 
from Eqs.(13) in terms of 9^ and 



as 
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V I (uj = E th + 0+)~V I [ — -1 
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~ s 2 M e ^(s 2 ) 9^, + s Mg^(s 2 ) gS\^\ = 5 g ^ 

-sMe^(s 2 )9 l p+M^(s 2 )gS\(j)\ = fy. (32) 

All Ms here can be obtained from those of the broad 
resonance together with the identification of |Ao| = <?|<^o| 
as follows: 



Me.es = M e&e& ; Mg^ = M 9&A + 2|</> 



M 



AA • 



(33) 



The behaviour of Ais at uj near the threshold energy can 
then be read off from the results in the broad resonance 
case above and the corresponding behaviour of the added 
terms. 

The presence of the constant inhomogeneous terms in 
Mg,$ and M^ in (33 1 allows for different damping be- 



haviour from that seen previously and we shall see that 
these differences are realized. 



A. The BCS regime 



For \i > in the BCS regime, together with Eqs.(17l, 
Mmr and Mg,<t,R follow 



M 



<t>4>R 



(w = E th + + )~Mt < t >R = 2(2 f i-i?); (34) 



with 



V R = M 
t> I =M 



, iM<p<pR . 



8M 2 



Ml 



(37) 



In particular, in addition to Mg^g^ /, T>j now becomes 
singular as u is near the threshold energy. Thus, the 
terms involving either Mg^g^i or T>i will become dom- 
inant in determining the late-time dynamics of the am- 
plitude mode. We then find that the integrands in the 
corresponding equation to Eq.(ll I, in terms of the vari- 
ables #0 and 6(f>, can be given by 



"DrMb^i -VlMg^g^R^ 



v 2 ~ 1 ' ' 



VrMs^i-VtM 

-,2 T T)2 \ui=E th +Q+ 



Mt 



;4>R I W 



(39) 



In general, the damping behaviour of the amplitude 
mode at late times then becomes 



x\mm M eM°) x 1 VRMe^i-ViMe^R 1 rp , ,. u 

9*\m - -^oT^-^ ^ ^^coe^t-*/^ 

+ 1 Mg^R _L ^ + 



That is, the decay of the amplitude mode in the narrow 
resonance follows the power-law i~ 3 / 2 , prior to reaching a 
saturated value in the BCS regime. This is different from 
the broad resonance result quoted earlier where the am- 
plitude mode decays dominantly as i^ 1 / 2 at late times in- 
stead. The only exception to this decay behaviour of the 
amplitude mode arises when the unitary limit (infinite 
scattering length) 2/i = v occurs in the regime of posi- 
tive chemical potential OHH]- We then have M.$$r(iS) 
of ( 33 1 vanishing linearly at u) near the threshold energy 

JSth, i.e. M^niuj = E th + 0+) w - lj, which gives 

the same £ -1 / 2 behaviour as its general counterpart in 



the broad resonance given in Eqs. (17). 

This finding of t -3 / 2 damping behaviour in the BCS 
regime seems to contradict the conclusion in [3] in which 
the same damping dynamics is found in both broad (one- 
channel model) and narrow (two-channel model) reso- 
nances. However, our idealization of narrow resonances 
for which we have taken the contact interaction strength 
U = in (fTl) cannot provide a valid description in the 
regime of deep BCS. For a more general two-channel 
model with a nonzero U, characterized by the back- 
ground scattering length, U e ff ~ U in the deep BSC 
regime [H [10] . Thus, the damping behaviour of [4] will be 
recovered in the deep BCS regime with a power law i -1 / 2 
decay, the same result as for the one-channel model, al- 
though the details are much messier. The outcome would 
be that, although this term might have a small coeffi- 
cient, it has the dominant long-time behaviour. It should 
be possible to see the otherwise subdominant <~ 3 / 2 be- 
haviour in the BSC regime if £/ «C g 2 /(2/i — v). 



B. The BEC regime 

In the case of negative chemical potential in the BEC 
regime, as a result of Eqs. (33) and the identification of 



|Ao | = sl^oli one finds that the corresponding Ms at w 
near the threshold energy bear similarity to those in the 
broad resonance in Eqs. (26), provided their coefficients 



are changed. The changes are summarized as follows: 



M^r(uj = E th + 0+) ~ M^r 
In. — TA 4- 

2tt 



2(2 M - D) + — 2 (2m M ) » ^ ^-L- 



Me^R^ = E th + 0+) ~ Me^R 
2^ 2V m) 2V2 ImI 



(41) 



(42) 



All other A4s follow Eqs. (26) and (27) by changing again 
the labels #a and A to 0a, and (j>. Additionally, |A | 



is replaced by <?|</>o|- Thus, XWx) follows Eq.(19) with 
T> R (j}, modified from Eqs. (28) and is thus given by 



V, 



9 \ c Po\ 



Then the d eca y of the amplitude mode at late times is 
given by Eq.(31), in which all Ads and Vs are replaced 



by those in terms of the variables 0^ and 6\<f>\, as follows: 



gs\<l>\(t) 



0(0) 



1 VbMb^i 



T>rM 6 



Vr 



1 V R Me^i-V I M e ^R 1 



E ^ tm 



cosfE'thi + tt/4](5 



cos^thi + 7r/4]<5e 



(44) 



The outcome, that the damping behaviour is t~ 3 ' 2 , as for 
broad resonances, is consistent with the argument in [3] 
but, yet again, there is the possibility of the condensate 
relaxing to a displaced value. 



V. CONCLUSIONS 

In this paper we determined the time evolution of a 
condensate of cold Fermi atoms, whose interaction can 
be tuned by a Feshbach resonance throughout the BCS 
and BEC regimes, in response to (small) non-equilibrium 
perturbations of the phase and amplitude modes. Inter- 
action between the collective modes and the constituent 



particles is key for our understanding of the relational 
dynamics, which occurs in a collisionless regime via the 
Landau damping effect. 

In particular, we have focussed on the evolution of the 
collective amplitude mode. For a broad resonance the 
amplitude mode decays as i -1 / 2 at late times in the BCS 
regime with positive chemical potential whereas it de- 
cays as i -3 / 2 in the BEC regime with negative chem- 
ical potential. Our conclusions, although obtained by 
different means, agree with those in the literature [3]. 
However, when the resonance is very narrow we disagree 
with [4], which sees no difference between the behav- 
ior for broad or narrow resonances. We find that for 
an idealized narrow Feshbach resonance the collective 
mode decays as i -3 / 2 throughout both the BCS and BEC 
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regimes. Nonetheless, we expect that the deviation from 
idealized narrow resonances restores the original longtime 
behaviour in the deep BCS regime. The i -3 / 2 behavior 
can possibly be seen experimentally in the BCS regime 
when the effect of the background scattering length can 
be largely ignored. More robustly beyond the results of 
[3] the condensate amplitude can decay to a value shifted 
from its initial condition, arising from the perturbation 
on the phase mode due to the amplitude-phase coupling. 
The results can be tested experimentally. 
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